Angular correlations in the two— photon decay of helium-like heavy ions 
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The two-photon decay of heavy, helium-like ions is investigated based on second-order pertur- 
' bation theory and Dirac's relativistic equation. Special attention has been paid to the angular 

^N) I emission of the two photons, i.e., how the angular correlation function depends on the shell struc- 

ture of the ions in their initial and final states. Moreover, the effects from the (electric and magnetic) 
non-dipole terms in the expansion of the electron-photon interaction are discussed. Detailed calcu- 
lations have been carried out for the two-photon decay of the 1825^5*0, ls2s'^Si and ls2p''Po states 
of helium-like xenon Xe^^"*", gold Au^^"*" and uranium U®°^ ions. 



PACS numbers: 31.10.-l-z, 31.30.jc, 32.80.Wr 



Ph; I. INTRODUCTION 

^ I Owing to the recent advances in heavy-ion accelerator and trap facilities as well as in detection techniques, new 
^ ■ possibilities arise to study the electronic structure of simple atomic systems in strong (nuclear) Coulomb fields. 
Relativistic, quantum electrodynamics (QED), or even parity non-conservation (PNC) effects, which are difficult to 
isolate in neutral atoms, often become enhanced in high-Z, few-electron ions. In order to improve our understanding 
of these fundamental interactions, a number of experiments have been recently carried out on the characteristic 
photon emission from heavy ions [l|-|3|. Apart from the one-photon bound-bound transitions, the two-photon decays 
of mctastable ionic states have also attracted much interest since the analysis of its properties may reveal unique 
QQ , information about the complete spectrum of the ion, including negative energy (positron) solutions of Dirac's equation. 
fNl ' Until now, however, most two-photon studies were focused on measuring the total and energy-differential decay rates 
which were found in good agreement with theoretical predictions, based on relativistic calculations [ll|-|20|. 
\^ . In contrast, much less attention has been previously paid to the angular and polarization correlations between the 
' emitted photons. The first two-photon correlation studies with heavy ions arc likely to be carried out at the GSI 
facility in Darmstadt, where the significant progress has been recently made in development of solid state, position 
sensitive x-ray detectors [l!], [2^. By means of these detectors, a detailed analysis of the angular and polarization 
properties of two-photon emission will become possible and will provide more insights into the electronic structure of 
, , ■ heavy, few-electron ions. 
r> ' . 

^ ' Despite of the recent interest in two-photon coincidence studies, not much theoretical work has been done so far to 
5^ I explore the photon-photon correlations in the decay of heavy atomic systems. While some predictions are available for 
the hydrogen-like [IS, [13] and neutral atoms , no systematic angular (and polarization) analysis was performed for 
the helium-like ions which are the most suitable candidates for two-photon investigations in high-Z domain. In the 
present work, therefore, we apply the second-order perturbation theory based on Dirac's equation to investigate the 
7 — 7 angular correlations in the decay of two-electron systems. The basic relations of such a (relativistic) perturbation 
approach will be summarized in the following Section. In particular, here we introduce the second-order transition 
amplitude that describes a bound-state transition under the simultaneous emission of two photons. The evaluation of 
these (many-body) amplitudes within the framework of the independent particle model (IPM) is thereafter discussed 
in Section fllBl Within this approximation that is particularly justified for the high-Z regime (25l - [27| . the photon- 
photon correlation function from Section fll Cl ean be traced back to the one-electron matrix elements. This reduction 
enables us to implement the well-established Green's function as well as finite basis set methods [28l - l30| and to 
calculate the correlation functions for the ls2s — > Is^ ^Sq, ls2s^Si — > Is^ and ls2p^Po — > Is^ ^S'o transitions 
in helium-like xenon Xe^^+, gold Au^^"'' and uranium U^^''' ions. Results from these computations are displayed in 
Section Hm and indicate a strong dependence of the photon emission pattern on the symmetry and parity of initial and 
final ionic states. Moreover, the significant effects that arise due the higher-multipole terms in the expansion of the 
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electron-photon interaction arc also discussed in the context of angular correlation studies. Finally, a brief summary 
is given in Section IIVI 



II. THEORETICAL BACKGROUND 



A. Second— order transition amplitude 



Since the second-order perturbation theory has been frequently applied in studying two-photon decay, here we 
may restrict ourselves to a short compilation of the basic formulas relevant for our analysis and refer for all further 
details to the literature [Til, [l^ H^. [la fisl. fiol [23 - [26l | . Within the relativistic framework, the second-order transition 
amplitude for the emission of two photons with wave vectors ki {i = 1, 2) and polarization vectors ux^ (A = ±1 ) is 
given by 



Mf.,{Mj,M,,X,,X-2) = 



IfJfMf 



7^■^(A;2,^tA2) 



IfJfMf 



7^t(A;2,■UAJ 



- E, + UJ2 



Ey — Ei + ui 



(1) 



where \jiJiMi) and \jfJfMf) denote the (many-electron) states with well-defined total angular momenta Jij and 
their projections Mij of the ions just before and after their decay, and jij all the additional quantum numbers as 
necessary for a unique specification. The energies of these states, Ei and Ef, are related to the energies wi_2 = cfci^2 
of the emitted photons by the energy conservation condition: 



Ei — Ef = hwi + tlhJ2 



(2) 



Using this relation, it is convenient to define the so-called energy sharing parameter y ~ wi/(a;i+cj2), i-C, the fraction 
of the energy which is carried away by the first photon. 

In Eq. ([1]), moreover, TZ is the transition operator that describes the interaction of the electrons with the electro- 
magnetic radiation. In velocity (Coulomb) gauge for the coupling of the radiation field this operator can be written 
as a sum of one-particle operators: 



TZ{k, Ux) = ^ dim Ax,m = ^ am ux e*'" '''" , 



(3) 



where ctm = ic(x,m,cty,m, cnz^m) dcnotcs the vector of the Dirac matrices for the m-th particle and .4a, m the vector 
potential of the radiation field. To further simplify the second-order transition amplitude ([T]) for practical compu- 
tations, it is convenient to decompose the potential Ax,m into spherical tensors, i.e. into its electric and magnetic 
multipole components. For the emission of the photon in the direction k = (9, cj)) with respect to the quantization 



(z-) axis such a decomposition reads [31 



uxe 



ihr 



^E E E *''W'^'(*^)'«'M(fc)^MA(fc): 
L=l M=-L p=0,l 



(4) 



where [L] = 2L + 1, D^j^ is the Wigner rotation matrix of rank L and a^^/ ' (k) refer to magnetic (p=0) and electric 
(p=l) multipoles, respectively. 

The multipole decomposition of the photon field in terms of its irreducible components with well-defined trans- 
formation properties enables us to simplify the second-order amplitude by employing the techniques from Racah's 
algebra. Inserting Eqs. ©-(jl]) into the matrix element ([T]) and by making use of the Wigner-Eckart theorem, we 
obtain: 

M/,(M/,M„Ai,A2) = 2^ J2 E (-0''^+^M^i,i2]^/'(-iAi)''M-*A2)^^^M\A,(fci)^ArA.(^2) 

LiMipi L2M2P2 
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E 



{JfMfLiAh\J,AU) {JM. L2M2\ J^M,) 5^V,i,p,(w2) 



+ {JfMfL2M2\JMu) {J.M,LiMi\J,Ah) 5^^^,, ^^^^ (c^i) 
where the second-order reduced transition amphtude is given by: 



(5) 



LlPl,L-2P2 



(^2) = E 



IfJf 



liJi 



En — Ei + UJ2 



(6) 



Here, the summation over the intermediate states formally runs over the complete spectrum of the ions, including a 
summation over the discrete part of the spectrum as well as the integration over the positive- and negative-energy 
continua. In practice, such a summation is not a simple task especially when performed over the many-electron states 
|7i/J^). In the next section, therefore, we shall employ the independent particle model in order to express the reduced 
matrix elements (jH]) for many-electron ions in terms of their one-electron analogs. 



B. Evaluation of the reduced transition amplitudes 

As seen from Eqs. ([5])-®, one has first to generate a complete set of many-electron states I7 J) in order to calculate 
the second order transition amplitude Mfi. A number of approximate methods, such as multi-configuration Dirac- 
Fock (MCDF) [s^,!!^ and configuration interaction (CI) [l5|, are known in atomic structure theory for constructing 
these states. Moreover, the systematic perturbative QED approach in combination with the CI method turned out to 
be most appropriate for describing both transition probabilities [13, HI] and transition energies [sl] in highly charged 
ions. In the high-Z domain, however, the structure of few-electron ions can already be reasonably well understood 
within the independent particle model (IPM). This model is well justified for heavy species especially, since the 
interelectronic effects scale with 1/Z and, hence, are much weaker than the electron-nucleus interaction (25l. [26l. ISTj . 
Within the IPM, that takes the Pauli principle into account, the many-electron wave functions are approximated 
by means of Slater determinants, built from one-particle orbitals. For this particular choice of the many-electron 
function, all the (first- and the second-order) matrix elements can be easily decomposed into the corresponding 
single-electron amplitudes. 

For a helium-like system, the decomposition of the reduced amplitude ([6]) reads: 



S 



J,. 

Llpi,L2P2 



jl' jo Jv 1 nj 



(7) 



where the one-electron matrix elements of the (electric and magnetic) multipole field operators are given by 



E 



(fc2 



El, — Ei + U!2 



(8) 



We assume here that the "spectator" electron, being in hydrogenic state |nojo), stays passive in the decay process. 
Moreover, \niji), \ni,j^) and \nfjf) = |noio) denote the initial, intermediate and final states of the "active" electron, 
correspondingly. The great advantage of formula ([T]) is that it helps us to immediately evaluate the many-electron 
transition amplitude ^ in terms of the (one-particle) functions Sf^^^^ j^^^^{u)2)- The summation over the complete 
one-particle spectrum that occurs in these functions can be performed by means of various methods. In the present 
work, we make use of (i) the relativistic Coulomb-Green's function j23 . [ssl . [s^ and (ii) a B-spline discrete basis set 
[20l [28l - l30| to evaluate all the second-order transition amplitudes. Indeed, both approaches yield almost identical 
results for the angular correlation functions in the two-photon decay of heavy helium-like ions. 



C. Differential decay rate 



Equation ([5]) displays the general form of the relativistic transition amplitude for the two-photon decay of many- 
electron ions. Such an amplitude represents the "building block" for studying various properties of the emitted 
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FIG. 1: (Color online) Angular correlation function pop for the ls2s^ Sq — > Is^ ^ So two-photon decay of helium-like xenon, 
gold and uranium ions. Calculations obtained within the electric dipole 2E1 approximation (dashed line) are compared with 
those including all the allowed multipoles (solid line). Results are presented for the relative photon energies y = Q.l (upper 
panel) and 0.5 (lower panel). 



radiation. For instance, the differential two-photon decay rate can be written in terms of (squared) transition 
amplitudes as: 



Aw UJ1UJ2 



dwidf^idfla (27r)3c3 2J, _ ^ ^ 

if we assume that the excited ions are initially unpolarized and that the spin states of the emitted photons remain 
unobserved in a particular measurement. As seen from expression ([9]), the two-photon rate is single differential- 
owing to the conservation law ^ — in the energy of one of the photons but double differential in the emission angles. 
Accordingly, its further evaluation requires to determine the geometry under which the photon emission is considered. 
Since no particular direction is preferred for the decay of an unpolarized (as well as unaligned) ion, it is convenient to 
adopt the quantization (2-) axis along the momentum of the "first" photon: fci||z. Such a choice of the quantization 
axis allows us to simplify the rate Q and to define the angular correlation function: 

W2, {e,y) = - Ef) ^^^J^^^^^ {01 = 0, = 0, h = 0), (10) 

that is characterized (apart from the relative energy y) by the single polar angle 6 — 62 oi the "second" photon 
momentum with respect to this axis. In this expression, moreover, the factor Stt^ arises from the integration over 
the solid angle dQ,i = sin^id^id^i of the first photon as well as the integration over the azimuthal angle d(/)2 of the 
second photon. In the next Section, we shall investigate the dependence of the function W2'y{0,y) on this opening 
angle 6 for various bound-bound transitions and for a range of (relative) photon energies. 



III. RESULTS AND DISCUSSION 



With the formalism developed above, we are ready now to analyze the angular correlations in the two-photon decay 
of helium-like heavy ions. In nowadays experiments, the excited states of these ions can be efficiently populated in 
relativistic ion-atom collisions. For example, the formation of the metastable \s2s^Sq state during the inner-shell 
impact ionization of (initially) lithium-like heavy ions has been studied recently at the GSI storage ring in Darmstadt 
[4(| . The radiative deexcitation of this state can proceed only via the two-photon transition ls2s ^5*0 — Xs^^Sq since 
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FIG. 2: (Color online) Angular correlation function pop for the ls2s'^Si — > Is'^ ^5*0 two-photon decay of helium-like xenon, 
gold and uranium ions. Calculations obtained within the electric dipole 2E1 approximation (dashed line) are compared with 
those including all of the allowed multipoles (solid line). Results are presented for the relative photon energies y = 0.1 (upper 
panel) and 0.5 (lower panel). 



a single-photon decay to the Is^ ^So ground state is strictly forbidden by the conservation of angular momentum. 
Fig. 1 displays the photon-photon angular correlation function for this experimentally easily accessible decay of 
helium-like xenon Xc^^+, gold Au^^+ and uranium 

U90+ 

ions and for the two energy sharing parameters y = 0.1 
(upper panel) and y = 0.5 (lower panel). Moreover, because the radiative transitions in high-Z ions are known to 
be affected by the higher terms of the electron-photon interaction ([3]), calculations were performed within both, the 
exact relativistic theory (solid line) to include all allowed multipole components {piLi, P2L2) in the amplitude ([5]) as 
well as the electric dipole approximation (dashed line), if only a single term with Li ~ L2 = I and pi ~ p2 = 1 is 
taken into account. In the dipole 2E1 approach, as expected, the angular distribution is well described by the formula 
1 + cos^ that predicts a symmetric — with respect to the opening angle 6 = 90° — emission pattern of two photons. 
Within the exact relativistic theory, in contrast, an asymmetric shift in the angular correlation function is obtained. 
As can be deduced from Eqs. (l7|)- (fT0|) . this shift arises from the interference between the leading 2E1 decay channel 
and higher multipole terms in the electron-photon interaction: 

W2^i9, y) oc (1 + cos2 e) + A^ cose-^^ cos^ + ... , (11) 

Sei 3 Sei 

where, for the sake of brevity, we have introduced the notation Slp — S'[p~l'p{LOi) + S'[p~l'p{uj2)- For high-Z domain, 
the photon emission occurs predominantly in the backward directions if the nondipole terms are taken into account; 
an effect which becomes more pronounced for the equal energy sharing (cf. bottom panel of Fig. 1). Including the 
higher multipoles into the photon-photon correlation function, a similar asymmetry was found in the past for the 
2s2/i ^ 1*1/2 decay in hydrogen-like heavy ions both within the nonrelativistic [2^ and relativistic p3 | theory. 

Apart from the singlet ls2s ^Sq, the formation of the triplet ls2s ^6*1 state has been also observed in recent collision 
experiments at the GSI storage ring [^,|4^. Although much weaker in intensity (owing to the dominant Ml transition), 
the two-photon decay of this ls2.s'^S'i state has attracted recent interest and might provide an important testing 
ground for symmetry violations of Bose particles (26l . |4l|. The angular correlation between the photons emitted 
in this ls2s'^5i — > Is'^^Sq (two-photon) decay is displayed in Fig. 2, by comparing again the results from the 
exact relativistic theory with the 2E1 dipole approximation. As seen from the figure, the photon-photon correlation 
functions for the 2'^ Si — > l^S'o transition is much more sensitive with regard to higher multipoles in the electron- 
photon interaction than obtained for the 2^5*0 l^S'o decay. The strongest non-dipole effect can be observed for 
the equal energy sharing {y — 0.5), where the two-photon emission is strictly forbidden within the electric dipole 
approximation. This suppression of the 2E1 decay is a direct consequence of the exchange symmetry of photons as 
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required by the Bose-Einstein statistics and, hence, aparticular case of the Landau- Yang theorem that forbids the 
decay of vector particles into two photons (cf. Refs. j26l l4l| - |43| for further details). In contrast to the 2E1 channel, 
the E1M2 2^Si — > I^Sq transition can proceed even if the energies of the two photons are equal. This transition as 
well as higher multipole terms give rise to a strongly anisotropic correlation function that vanishes for the parallel 
and back-to-back photon emission and has a maximum at 9 = 90°. 

Large effects due to the higher multipole contributions to the ls2s'^Si Is'^-^Sq two-photon transition can be 
observed not only for the case of equal energy sharing (y ~ 0.5). For the relative energy y ~ 0.1, for example, 
the photon-photon angular correlation function is found symmetric with regard to 6 = 90° in the electric dipole 
(2E1) approximation but becomes asymmetric in an exact relativistic theory. In contrast to the decay of 2^5*0 state, 
however, a predominant parallel emission of both photons appears to be more likely if the higher multipoles are taken 
into account. For the 2^51 — >■ 1^5*0 two-photon decay of helium-like uranium U^°+, for example, the intensity ratio 
W2'y{9 = 0°, y = O.1)/W2-y(0 = 180°, y = 0.1) increases from unity within the electric dipole approximation to almost 
1.6 in the exact relativistic treatment. 

Until now we have discussed the photon-photon correlations in the decay of ls2s (singlet and triplet) helium-like 
states. Besides these — well studied — transitions, recent theoretical interest has been focused also on the ls2p^Po — 
Is^ ^5*0 two-photon decay whose properties are expected to be sensitive to (parity violating) PNC phenomena in atomic 
systems (27| . Future investigations on such subtle parity non-conservation effects will require first detailed knowledge 
on the angle (and polarization) properties of two-photon emission as well as the role of non-dipole contributions. 
The angular correlation function ([TU]) for the 2'^Po ^ l^S'o transition is displayed in Fig. 3, again, for two relative 
photon energies y = 0.1 and 0.5 and for the nuclear charges Z = 54, 79 and 92. Calculations have been performed 
both within the exact theory and the ( "electric and magnetic" ) dipole approximation which accounts for the leading 
ElMl-MlEl decay channel. As seen from the figure, the emission pattern strongly depends on the energy sharing 
between the photons. If, for example, one of the photons is more energetic than the second one their parallel emission 
becomes dominant (cf. upper panel of Fig. 3). In contrast, photons with equal energies, i.e. when y = 0.5, are more 
likely to be emitted back-to-back while the differential rate © vanishes identically for 9 = 0°. Such a behaviour of 
the photon-photon angular correlation function is caused by the interference between two pathways which appear for 
each multipole component of the 2^Pq — > 1^5o transition. For instance, the leading ElMl-MlEl decay may proceed 
either via intermediate n^Si or n^Pi states, thus given rise to a "double-slit" picture that becomes most pronounced 
for the equal energy sharing. Simple analytical expression for the angular correlation function which accounts for 
such a Young-type interference effect can be obtained from Eqs. ©-(O as: 



W2y{9,y) cx sm''9/2\SEiMif f 1 + 2(1 + 2 cos 0)^^HMi^ 



+ cos4 0/2 \Veimi\^ (l-2{l-2cos9)^^^\ + (12) 

V J-^EIMI J 

where, similar as before, we denote ^Lpi.Lp^ = Sl"=\p2(^i) + •^'i^'ips ('^^2) + ^'l^^, ipi ('^z) + -S'i^^^ip^ (wi) and 
"DlpuLp, = Sipl'ip^ (uJi) - Sip=\j,^ (W2) + (^2) - ^^;r.ipi (^i)- Obviously, if the energies of the two photons 

are equal, wi = W2, the second term in Eq. (fT2|) turns to be zero and the photon emission is described by the sin^ 9/2 
angular distribution modified by the non-dipole terms in the expansion of electron-photon interaction. As seen from 
the lower panel of Fig. 3, the contribution from these terms becomes more pronounced for the back-to-back photon 
emission (6* = 180°) where they lead to about a 30 % enhancement of the correlation function. It is interesting to 
note that such an enhancement remains almost constant along the helium isoelectronic sequence for Z > 54 due to 
similar (cx Z^^) scaling of the ElMl and E2M2 transition probabilities. Therefore, our calculations clearly indicate 
the importance of higher multipoles for analyzing the photon-photon correlations not only for high-Z domain but 
also for medium-Z ions. 



IV. SUMMARY AND OUTLOOK 



In summary, the two-photon decay of heavy, helium-like ions has been investigated within the framework of the 
relativistic second-order perturbation theory and the independent particle model. In this study, special emphasis 
was placed on the angular correlations between the emitted photons. A general expression for the photon-photon 
correlation function was derived that accounts for the complete expansion of the radiation field in terms of its multipole 
components. Based on solutions of Dirac's equation, this function has been calculated for the two-photon decay of the 
ls2s ^^o, ls2s ^5*1 and 1s2p^Pq states of helium-like xenon Xe^^+, gold Au^^+ and uranium 11^"+ ions. As seen from 
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FIG. 3: (Color online) Angular correlation function pop for the ls2p ''Po Is^ ^Sq two-photon decay of helium-like xenon, 
gold and uranium ions. Calculations obtained within the dipole ElMl approximation (dashed line) are compared with those 
including all of the allowed multipoles (solid line). Results are presented for the relative photon energies y = 0.1 (upper panel) 
and 0.5 (lower panel). 



the results obtained, the photon emission pattern appears to be sensitive to the symmetry and parity of the particular 
excited state as well as to the higher multipole contributions to the electron-photon interaction. The strongest non- 
dipole effects have been identified for the ls2s^Si — >■ 1 ^Sq two-photon transition for which the 2E1 decay channel 
is forbidden owing to symmetrization properties of the system. For the other two transitions, ls2s ^5*0 — >■ 1 and 
ls2p^Po 1 ^Sq, the higher multipoles of the radiation field typically result in a 10-30 % deviation of the photon- 
photon correlation function from the (analytical) predictions obtained within the dipole 2E1 approximation. This 
deviation becomes most apparent for the parallel and back-to-back photon emission and may be observed not only 
for high-Z but also for medium-Z ions. 

The second-order perturbation approach based on the independent particle model, used in the present calculations, 
is appropriate for the analysis of forthcoming experimental studies on the two-photon transitions between the ^'^'^^Lj 
excited and the ground states of helium-like, heavy ions. Besides these spontaneous decays, whose energies usually 
reach 100 keV, induced J = 0— >-J = + 27 transitions between excited states are also likely to be explored at the 
GSI ion storage ring Having energies in the optical range (2-3 eV), these transitions may provide an alternative 
and very promising tool for studying the parity violation phenomena. Their theoretical analysis, however, requires a 
more systematic treatment of the electron-electron interaction effects. Based on the multi-configuration Dirac-Fock 
approach and B-spline basis set method, investigations along this line are currently underway and will be reported 
elsewhere. 
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